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Abstract 

We consider the Hall conductivity of two-dimensional non-interacting Bloch 
electrons when the magnetic flux per unit cell is a rational number p/q where 
p and q are mutually coprime. We present a counter-example for the naive 
expectation that the Hall conductivity carried by a band is given by treating 
gap minima as Dirac fermions. Instead of the above expectation, we show that 
the change of the Hall conductivity at a gap-closing phenomenon is given by 
the Dirac fermion argument. Comparing with the Diophantine equation, our 
result implies that a band-gap closes at q points simultaneously. Furthermore, 
we show that the dispersion relation is (/-fold degenerate in the magnetic 
Brillouin zone. 

PACS numbers: 73.40.H, 03.70, 02.40.P 



Typeset using REVTgX 



*To appear in Phys. Rev. B 

^Present Address: Department of Applied Physics, University of Tokyo. Hongo, Bunkyo-ku, 
Tokyo 113 JAPAN. 

* E-mail address: oshikawa@mmm.t.u-tokyo. ac.jp 



I. INTRODUCTION 



It was a surprise that Hall conductance of two-dimensional electron system was found 
to be quantized in some experiment s.0 The presence of disorder or edges are considered to 
be necessary for the quantization of the Hall conductance! in the experiments. However, 
basic important facts are that there are energy gaps in the bulk system without disorder 
and that the Hall conductivity is quantized when the Fermi level lies in the gap. In the free 
electron case, it is easy to show this; eigenstates are the degenerate Landau levels separated 
by energy gaps, and each Landau level contributes e 2 /h exactly to the Hall conductivity. 

When a periodic potential is present, it is a more non-trivial problem Since a Landau 
level splits into several subbands by turning on a periodic potential, one may expect that a 
subband carries fractional (in units of e 2 /h) Hall conductivity. However, Thouless, Kohmoto, 
Nightingale and den Nijsi (TKNN) showed that the Hall conductivity carried by a magnetic 
subband is always an integer. This quantization comes from the topological nature; TKNN 
integer is a topological invariant on the magnetic Brillouin zone!. Ishikawa et al% also 
discussed the topological aspect of the quantized Hall conductivity. Topological character 
of the Hall conductivity carried by edge states and its relation to the bulk TKNN integer 
are also discussed!. 

On the other hand, it has been widely accepted that low-energy behavior of a thermo- 
dynamic system is correctly described by a continuum field theory. In the present case, a 
minimum of a gap between magnetic bands may be described by a Dirac fermion (in 2+1 
dimensions), which gives the Hall conductivity — e 2 /2h sgnm ( m is the Dirac mass.) Hence 
we may expect that the Hall conductivity carried by a band is given by counting only the 
Dirac fermions. There are several wor based on this idea. 

However, it is not trivial whether the Dirac fermion argument is true or not, because the 
Hall conductivity is given by an integral over whole Brillouin zone while the Dirac fermion 
argument looks only the gap minima in the Brillouin zone. Therefore, in this paper we 
investigate the validity of the Dirac fermion argument for Bloch electrons, and its relation 
to the topological TKNN integer. 

Contrary to the naive expectation, we found that the Hall conductivity carried by a band 
cannotgiven by the Dirac fermion argument in general. Instead, applying the argument by 
SimonO and otherll, we showed that the change of the Hall conductivity at a gap-closing 
point is correctly described by the Dirac fermion argument. 

Comparing with the Diophantine equation, this result implies that a band-gap closes 
simultaneously at q points. Moreover, we showed a stronger statement that the dispersion 
relation is g-fold degenerate in the magnetic Brillouin zone for general 2D Bloch electrons. 
This is a generalization of the result0 for the tight-binding model on a square lattice. 

Throughout this paper, we consider two-dimensional non-interacting electron systems 
under a periodic potential and a uniform magnetic field at zero temperature. The flux per 
unit cell is assumed to be a rational number, and the Fermi level lies at a band gap except 
when the gap closes. We set the velocity of light c to be 1 while we keep Planck's constant 
h or h = h/2n. 

The organization of this paper is as follows. In Section [H], we give a brief review on the 
Hall conductivity as a topological invariant. In Section |TJ, we review the Dirac fermion in 
2 + 1 dimensions. Section fV| describes a counterexample in which Dirac fermion argument 
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cannot give the correct value of the Hall conductivity of Bloch electrons. In Section |V|, we 
show that the change of the Hall conductivity during a gap-closing is given by the Dirac 
fermion argument. A g-fold structure in the magnetic Brillouin zone is shown in Section [VT|. 



Conclusions and a discussion are given in Section |VII . 



II. HALL CONDUCTIVITY OF BLOCH ELECTRONS 

Here we review known important facts about the quantized Hall conductivity of Bloch 
electrons, especially the topological aspectiS of the quantized Hall conductivity. 

The Schrodinger equation for a two-dimensional non-interacting electron system in a 
uniform magnetic field is written as 

Hij){f) = [^-(P + eA) 2 + U{r)]ip(r) = E^j(r) (1) 

where p is the momentum —iTiW and A is the vector potential. In our case, the potential 
U(r) is periodic, i.e., 

U(f + aO = U(f + a 2 ) = U(r) (2) 

where a\ and a 2 are linearly independent Braveis vectors. We take the symmetric gauge 
A = 1/2(B x r) for the moment. Let us define the magnetic translation operators 

f li = exp[ t -R-(p-eA)]. (3) 

In the symmetric gauge, this operator commutes with the kinetic term in the Hamiltonian, 
as well as with the potential term. 

We consider the case where the magnetic flux per unit cell is p/q, where p and q are 
mutually prime integers. Namely, the magnetic flux density B satisfies 

S.(^xa 2 ) = --. (4) 
q e 

Then the following relations hold: 

[T gSl ,T S2 ]=0,[T R ,H} = (5) 

where R = n(gai) + ma^ and n,m are integers. Thus we can apply the Bloch's theorem if 
we take an enlarged unit cell (magnetic unit cell) which is q times larger than the original 
unit cell. Correspondingly, the reciprocal space (magnetic Brillouin zone) becomes 1/q of 
the original Brillouin zone. Namely, the Schrodinger equation can be reduced as 

H{k)ui{r) = ^4(r) (6) 
H(k) = —(-ihV + hk + eA) 2 + U{r) (7) 
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where i is the band index, k is the crystal momentum and u l is a wavefunction with the 
generalized Bloch condition: 

u l {f) = u l {f+ qcii) exp [— ^l^B ' (^* x ^1)] 

6 -» 

= w*(f +a 2 )exp[-i— 5- (f x a 2 )]. (8) 

Defining the reciprocal lattice vectors g t by ^ • Oj = 5^, we have equivalence relations 

2tx 

k~k + —{fx ~ k + 2tT(?2 (9) 

— * — * 

among crystal momenta fc. Hence we can restrict k to the magnetic Brillouin zone 

-> 7T 71 

k = fcl<7l + ^2^2 ( < < -, "TT < k 2 < 7r). (10) 

From the linear response theory (Nakano-Kubo formula), the Hall conductivity of the 
system is given by 

v e 2 1 



i\Ei<Ep h 2m 



X 



J d k Lvc [{ ~dk~ r ~dkV ~ { ~dk~ r ~dhT ] - 



That is, the total Hall conductivity a xy is obtained as a sum of the contributions from all 
bands below the Fermi level. 

The contribution from a single band can be written in a compact form as follows. We 
will omit the band index when a single band is considered. For each band, we define a vector 
field in the magnetic Brillouin zone by 

A(k) = (u(k)\V k \u(k)) = J Muc d 2 rul{r)V k u % {r) (12) 

where is a vector operator (d/dk x , d/dk y ) and MUC represents the magnetic unit cell. 
It should be noted that we can choose the phase of the wavefunction for each k arbitrarily. 
Transformation of the phase as 

u' % (f) = u n {f)e^\ A\k) = A{k) + iV k f{k) (13) 

is nothing but a U(l) gauge transformation; A(k) is a U(l) gauge field defined on the 
magnetic Brillouin zone. Using this gauge field, the Hall conductivity carried by a single 
band is written as 



e 2 



1 f d ^kV k x A(k) (14) 
h 2m Jmbz 



where V&X denotes the rotation in two-dimensional /c-space and MBZ represents the mag- 
netic Brillouin zone. 
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Although a naive application of the Stokes' theorem implies that (|14]) is zero, a non- 
zero value of the hall conductivity arises from a non-trivial topology of the fiber bundle. 
That is, one cannot determine the phase of the wavefunction smoothly and uniquely 
over the entire Brillouin zone in general. For example, one can fix the gauge by requiring 
that the amplitude (a\u(k)) is real (|a) is an arbitrary wavefunction.) Though this seems 
a well-defined gauge fixing, this cannot fix the gauge at the zeros of the amplitude. We 
must cover the region near a zero by another patch, in which a different gauge fixing is 
chosen. Phase mismatch between the two gauges produces a non-zero value of (|T4]). We call 
a zero of the amplitude as a vortex. Uniqueness of the phase in each patch implies that the 
phase mismatch around a vortex should be an integral multiple of 2tt. This integer is called 
vorticity. Thus the Hall conductivity carried by a band ( |14|) is e 2 /h times the total vorticity 
in the magnetic Brillouin zone. While the location of the vortices depends on the gauge, 
total vorticity is gauge invariant and is known as the first Chern number of the principal 
17(1) bundle. 

III. DIRAC FERMION IN 2+1 DIMENSIONS 

To fix the convention, here we briefly review the relativistic Dirac fermion in 2 + 1 
dimensions. The Lagrangian density is given by 

C = i){ihd^ - eA^ip - m^) (15) 

Where is the (background) vector potential. The 7-matrices satisfy the anticommutation 
relation {7^,7^} = rf v where if" 1 " is the flat Minkowski metric. In the 2+1 dimensions, 7- 
matrices are 2x2 matrices. We choose the convention 

7 ° = o-^y = ta y n 2 = -%a x (16) 

where a x,y ' z are the standard Pauli matrices. When the electromagnetic field is absent, the 
1-body Hamiltonian is derived from (O) as 

H = ma z + p x a x + p y a v (17) 

where (p x ,P y ) represents the momentum. 

In the low-energy limit, the current is evaluated by several methods01ll 

(ef) = e^r^P) = —tr'FnSgn m (18) 

where F va is the field strength d v A a — d a A v . Namely, the Hall conductivity carried by a 
Dirac fermion is — e 2 /2/isgn m. 



IV. FAILURE OF THE DIRAC FERMION ARGUMENT 

In many cases, the low-energy behavior of a thermodynamic system can be treated by 
a continuous field theory. From the field-theory point of view, this means a field theory is 
independent of details in the regularization. 
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In the case of the Bloch electron system in a magnetic field, the low-energy states near 
a gap minimum may be treated as a Dirac fermionlii Thus we may expect that summing 
up contributions from each Dirac fermion gives the Hall conductivity of a band. Several 
authors made discussions based on this idea, and in fact it seems correct in some simple 
lattice mo delJl. 

However, we show that this naive expectation is false in general. As a counterexample, we 
take a tight-binding model discussed in Ref. [TJ]. We consider a isotropic square lattice with 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) hoppings with 1/3 magnetic flux 
per unit cell. Let the absolute value of the NN and NNN hoppings be 1 and t c , respectively. 
Namely, the Hamiltonian of the system reads 

H =Y1 C l+l,mCn,m exp (t6 A ) + C^m+l^m exp (i6 B ) 
n,m n,m 

-\-tc C n+l,m+l C n,m exp (i0c) 
n,m 

(i6 D ) + H.c. 

n,m 

where 

2-7T 2tt 1 

0a = O, 6 B = —n, c = D = —(n + -). (19) 

The energy spectrum of this model is analyzed by Hatsugai and KohmotoB. According to 
them, the energy eigenvalue E for each (crystal) momentum k is determined by the equation 

F(E) = f(k) (20) 

where F(E) is some polynomial of E and f(k) is given by 

f(k) = 2 cos (3^) [1 - 3t c 2 ] + 2 cos (3k y ) [1 - 3t 2 ] 

-2t c 3 {cos [3{k x + k y )} + cos [3{k x - k y )]}. (21) 

Precisely speaking, a minimum of a energy gap cannot always be regarded exactly as a 
Dirac fermion. In general Bloch electron systems, there are several bands separated by the 
energy gaps. If the gap is very small compared with other gaps, we can ignore other bands. 
Then the low-energy states near a gap minimum can be considered as one-particle states 
of a Dirac fermion. It is possible that a gap-closing point become singular; in this case the 
dispersion of the low-energy states is not relativistic. We have checked this is not the case in 
our counter-example. We will make a detailed discussion on this point in the next section. 

In our model, the first gap closes simultaneously at three points when t c ~ 0.268. An 
extremal point of an energy band is given by an extremal point of f(k), and it is easy to show 
that the first gap has three minima in the neighborhood of the gap-closing point t c ~ 0.268. 
Thus there are just three Dirac fermions in this region, if the Fermi level lies in the first 
gap. The dispersion relation of the model t c = 0.25, which is near the gap-closing point, is 
shown in Fig. 1. 

Since a Dirac fermion contributes ±e 2 /2h to the Hall conductivity a xy , it should be half- 
odd-integral multiple of e 2 /h from the Dirac fermion argument. However it should be an 
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integral multiple of e 2 /h according to the general theory of TKNN. Hence the naive Dirac 
fermion argument is incorrect. In other words, there is no 'anomaly-cancelling partner'! in 
our model. 



V. CHANGE OF THE HALL CONDUCTIVITY AND THE DIRAC FERMION 

The naive expectation that the Dirac fermion argument gives the Hall conductivity 
carried by a band is disproved in the last section. Nevertheless, the Dirac fermion argument 
still makes sense for our problem. We proved the following proposition. 

When a parameter in the Hamiltonian is varied, the Hall conductivity changes 
only where the gap closes. This change of the Hall conductivity is correctly given 
by the Dirac fermion argument; the Hall conductivity decreases by e 2 /h if the 
Dirac mass changes from negative to positive, and increases by e 2 /h if it changes 
from positive to negative. If the gap closes at several points simultaneously, the 
change of the Hall conductivity is given by summing up the contributions from 
each gap-closing point. 

The above proposition is shown by mapping of the low-energy states to the Dirac fermion 
and calculation of the change of the Chern number. The important point is that the change 
of the Chern number is determined only by the neighborhood of the gap-closing point. 
This fact was noticed by SimonS, and also by Avron et ai.@ in the context of the network 
problem. We also note that the proof has similar structure to the "intuitive topological 
proof "0 of the Nielsen-Ninomiya theorem in 3 + 1 dimensions. 

To discuss the change of the Hall conductivity, the Hamiltonian is assumed to vary 
smoothly. Namely, we discuss a family of Hamiltonians labelled by a parameter A. The 
reduced Hamiltonian H also should be parameterized by A in the reduced Schrodinger equa- 
tion (H). Thus H depends on three parameters, k x , k y and A; it can be said that H is defined 
in a three-dimensional space whose coordinates are k x , k y and A. 

First let us focus on the gap-closing phenomena. Assume that there is a gap-closing point 
(k*,X*) where energies of two bands are degenerate. Here we consider the generic case in 
which only two bands touch at this point. We can neglect other bands near this gap-closing 
point. Furthermore, we expand the Hamiltonian to the first order in 

P = (Vx-.Vy.Vz) = (h(k x - k x *), h(k y - k y *), A - A*) (22) 

where p is a three-dimensional vector. Since the effective Hamiltonian for the two bands is 
2x2 Hermitian matrix for each p, the most general form of the expansion is 

H = E* + b ■ p + V»o» V v + 0(p 2 ) (23) 

where /i, v = 1,2,3 and a 1,2 ' 3 = a x ' y,z . If V is singular (det V = 0), the dispersion relation 
near the gap-closing point is not Dirac-fermion like. We don't consider such non-genericll 
cases in this paper. 

In order to see this is a Hamiltonian of the Dirac fermion, we make a unitary transfor- 
mation and a redefinition of (v x ,Vy)- Since the term proportional to the identity matrix 
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can be absorbed to the redefinition of the energy and hence irrelevant, we only consider the 
terms containing a Pauli matrix. It is noted that a SU(2) transformation U induces the 
transformation on the Pauli matrices as 

Ua" U- 1 = W v o v (24) 

where R% is a 5*0(3) matrix. There is always an appropriate SU{2) matrix U corresponds to 
some given R. Any regular matrix V can be decomposed as (Gram-Schmidt decomposition) 

V = OT (25) 

where O is an orthogonal matrix and T is an upper triangle matrix with positive diagonal 
elements. Let us choose = OW where W = diag(l, 1, sgn(det V)). Then U transforms 
the Hamiltonian as 

UHU- 1 ~ o»TPW v p p v (26) 

where the terms proportional to the identity matrix are omitted. Let us define 

V» = (PxiPyiPz) = T u ^p u . (27) 

Since T is an upper triangle matrix with positive diagonal elements, this is a parity- 
conserving Affine transformation on p x ,p y and a scale transformation on p z : p z = 
(positive constant) x p z . 

By this redefinition, the Hamiltonian becomes 

H ~ sgn(det V)p z a z + p x a x + p y a y , (28) 

which is nothing but the Dirac Hamiltonian flT7| ) with the mass sgn(det V)^. 

In this way, the neighborhood of a gap-closing point can be regarded as a Dirac fermion 



and the sign of the mass is opposite before and after the gap-closing. In Ref. |17|, the 
Hamiltonian which is formally same as ([58]) represents the one-body Weyl Hamiltonian in 
3+1 dimensions. In this case, p z denotes the third momentum, and the Weyl fermion is left- 
handed (right-handed) if det V > (det V < 0). The similarity between our argument and 
that of Ref. |17] is based on the formal similarity between the one-body Dirac Hamiltonian 
in 2+1 dimensions and the one-body Weyl Hamiltonian in 3+1 dimensions. However, for 
our proof we should be careful not to mix the parameter A with the momenta k x , k y during 
the transformation. This point has been solved by the Gram-Schmidt decomposition. 

However, in our problem, we should be careful not to mix the parameter A with the 
momenta k x , k y during the transformation because each section for constant A represents a 
physical model. This point has been cleared by use of the Gram-Schmidt decomposition (pop. 

Next, we calculate the change of the Chern number as discussed in Refs. p|[12| . We 
introduce vortex lines in the three-dimensional parameter space. The vortex lines are defined 
by 

{(k x ,k y ,X)\3j:(a\u^k;X)) = 0} (29) 
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with respect to some reference vector (in the reduced Hilbert space) \a). The real part and 
the imaginary part of the above condition fix two degrees of freedom. Thus (^9|) defines a 
set of curves (vortex lines) in the three-dimensional parameter space. When a vortex line 
satisfies (a\u l (k; A)) = for the band i, the vortex line is said to be in the z-th band. 

It should be noted that the definition of the vortex lines is similar to that of the vortices. 
In fact, a section of the three-dimensional parameter space for a certain constant A represents 
the magnetic Brillouin zone, and its intersection with the vortex lines are vortices. We define 
the orientation of a vortex line as shown in Fig. 2. it is taken in such a way that A increases 
(decreases) if the vorticity at the section is positive (negative). 

The following propositions^ are central for our result. 

• For each gap-closing point, there is always a vortex line passes through it. 

• The vortex line passes through the gap-closing point upward (i.e. from the lower band 
to the upper band) along the defined orientation if det V > and downward if det V < 
0. 

Now we are going to examine the change of the Hall conductivity at the gap-closing point. 
We assume that the Fermi energy always lies in the gap except when it closes. We consider 
the case det V > 0. Let us consider the vortex line which passes through the gap-closing 
point. If the orientation of the vortex line at the gap-closing point is in the direction that A 
increases, the local vorticity is +1. According to the above proposition, the vortex line lies 
in the lower band when A < A* and in the upper when A > A*. (See Fig. 3.) 

If the orientation of the vortex line at the gap-closing point is the reverse, i.e. in the 
direction that A decreases, the vorticity is —1. Thus in this case the vortex line lies in the 
upper band when A < A* and in the lower when A > A*. When A is increased through the 
gap-closing value A*, the total vorticity of the lower band decreases by 1 in the cases; the 
Hall conductivity decreases by e 2 /h irrespective of the direction of the vortex line. 

On the other hand, the mass of the Dirac fermion is negative when A < A*. If A is 
increased, the mass become positive when A > A*. The prediction of the Dirac fermion 
argument is that the Hall conductivity decreases by e 2 /h when A is increased through the 
gap-closing value A*. Now we can see it gives the correct change of the Hall conductivity. 
Similar arguments can be applied for a gap-closing point with det V < 0. It is also easy to 
see that the change of the Hall conductivity is given by summing up the contributions from 
each gap-closing point when the two bands touch simultaneously at several points. 

We also note that, the vortex only moves to another band and never appear or disappear 
when two bands touch. This implies that the sum of the Hall conductivity carried by two 
bands is conserved in a gap-closing. This is a simple illustration for the conservation law 
first discussed by Avron, Seiler and SimonEl. 

Let us show how our result applies to the example presented in Section [TV[ The first gap 
closes simultaneously at three points when t c = t* ~ 0.268. The matrix V for these points 
can be calculated numerically and we obtain det V ~ 7.2 > for each point (see also the 
next section). Thus the Dirac mass is negative when t c < t* and positive when t c > t*, and 
the Hall conductivity of the fist band decreases by 3e 2 /h when t c is increased through t*. 
This is consistent with the result a xy changes from e 2 /h to —2e 2 /h obtained by Hatsugai 
and Kohmoto@. 
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VI. PERIODIC STRUCTURE IN THE MAGNETIC BRILLOUIN ZONE 



It is known that the Hall conductivity of the Bloch electrons satisfies the Diophantine 
equationl'El 

r = qs r + pt r (30) 

where the Fermi level lies at the r-th gap, s r is an integer and the total Hall conductivity 
is t r e 2 /h. A gap-closing does not affect the value of p,q and r. Thus if t' r e 2 /h is the Hall 
conductivity after the gap-closing, we have 

q(s r -s' r )+p(t r -t' r ) = (31) 

for an integer s' r ; the change of the Hall conductivity should be an integral multiple of qe 2 /h. 

Comparing with this, our result implies that two bands touch simultaneously at multiple 
of q points. In fact, in our example shown in Section [V the first gap closes simultaneously 
at q = 3 points. For the tight-binding model on a square lattice, KohmotcEl showed that 
the dispersion relation is g-fold in the magnetic Brillouin zone, combining a duality trans- 
formation^ and a gauge transformation. Although it seems difficult to extend the duality 
transformation to general models, we show that the g-fold structure is common in general 
Bloch electron systems. 

For convenience, we make a gauge transformation in the reduced Schrodinger equation (^|) 
to a Landau-like gauge as 

p 

u'g(r) = exp [-in-(f- <fi)(f ■ gfijutff) (32) 

A' = -{B x f) - ~[f ■ (B x a 2 )g 2 — r ■ (B x a^] 

= [{B x f) • a 2 \g 2 . (33) 

If we shift the momentum by k — > k + 2n(p/q)g 2 , the kinetic term changes as 

— f-iW + hk + eA) 2 -> — i-ihV + hk + h-g 2 + eA) 2 . (34) 
2m 2m q 

This can be absorbed by the shift r — > f — Ex, which changes the vector potential as 

A -> A - [(ai x a 2 ) • B]g 2 

= A-- V -g 2 (35) 
e q 

This shift does not change the potential term. Thus we see that E{k + 2ii( y p/q)g 2 ) = E{k) 
Since p and q are coprimes, for an arbitrary integer m there is always an integer n so that 

np = m (mod q). (36) 

This leads to the symmetry 
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_,,r* Z7T77J _ . 

£(£; + g 2 ) = E{k). 



(37) 



Hence the above symmetry means that the dispersion relation is g-fold; the magnetic Bril- 
louin zone fllOl ) consists of q sub- zones which have the same dispersion relation. 

The fact that two bands should touch simultaneously at multiple of q points follows from 
this proposition. Moreover, we can also see that the Dirac masses (if they are well-defined) 
are the same for the q Dirac fermions connected by the shift k — > k + 2n(m/q)g2- 

The above result shows the consistency between the Dirac fermion argument and the 
Diophantine equation. We note that, in general, the Diophantine equation ( p0[ ) implies that 
the Hall conductivity is not an integral multiple of qe 2 /h nor qe 2 /(2h), although the disper- 
sion is g-fold. Hence a Dirac fermion argument cannot give the correct Hall conductivity 
as long as one defines Dirac fermions with respect to the dispersion relation. The g-fold 
structure restricts only the change of the Hall conductivity to be a multiple of qe 2 /h. While 
the Dirac fermion argument can determine the precise value of the change of the Hall con- 
ductivity, the Diophantine equation restricts also the value of the total Hall conductivity; 
they are consistent and complementary results. 

VII. CONCLUSIONS AND DISCUSSION 

We discussed two-dimensional non-interacting Bloch electrons in a uniform magnetic 
field. The naive expectation that the Hall conductivity carried by a band is given by treating 
gap minima as Dirac fermions was found to be false. Instead of the naive expectation, we 
showed that the change of the Hall conductivity when the band gap closes is correctly 
given by the Dirac fermion argument. Comparing with the Diophantine equation, our result 
implies that the gap-closing occurs simultaneously at multiple of q points. We proved a 
stronger statement that the magnetic Brillouin zone consists of q sub-zones that have the 
same dispersion relation. 

Although there is no reason that a naive field-theory prediction should be always true, one 
may ask why it fails in the present case. Our answer is that the failure is already implicit 
in the (naive) Dirac fermion argument itself. The Hall conductivity of a Dirac fermion 
depends only on the sign of the Dirac mass. That is, a Dirac fermion with arbitrary large 
mass contributes to the Hall conductivity. However, the large Dirac mass corresponds to the 
large gap between the bands in the Bloch electrons. When the gap is not small compared 
with other gaps, the identification between the Dirac fermion and the gap minimum becomes 
ambiguous. 

On the other hand, since the Hall conductivity is given by an integral over the whole 
magnetic Brillouin zone, it seems impossible to determine the value only from gap minima. 
However, the Hall conductivity is not an ordinary integral but a topological invariant and 
cannot change except when the gap closes. Hence we can expect that the change of the 
Hall conductivity can be described by the local neighborhood of the gap-closing point. This 
expectation is realized in the proof. 

Finally, we comment on the Dirac fermion argumenti using the Widom-Streda formula. 
It is summarized as follows. The Hall conductivity is related to the change of the charge 
density in an infinitesimal extra magnetic field. In a constant magnetic field, the energy 
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levels of a Dirac fermion form Landau levels. There are zero modes, whose energy depends 
on the Dirac mass, among the Landau levels. Since the spectrum besides the zero modes 
is symmetric for the Dirac fermion, the Hall conductivity is determined by the zero modes 
as — e 2 /(2/i)sgnm for each Dirac fermion. However, in a general Bloch electron system, the 
spectrum is not symmetric about the Fermi level; the Dirac fermion argument is not reliable 
to obtain the Hall conductivity. Nevertheless, we can expect that only the contribution from 
the zero modes will change at a gap-closing point. This is consistent with our result that 
the Dirac fermion argument gives only the correct change of the Hall conductivity. 
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Figure 1: The energy bands for the whole magnetic Brillouin zone of the 
tight-binding model on the square lattice with the NN and NNN hoppings. 
Here the NNN hopping t c — 0.25, which is slightly less than the gap-closing 
point t c ~ 0.268. If the Fermi level lies in the first gap, there are three 
(approximate) Dirac fermions corresponding to the minima of the gap. 
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Figure 2: The definition of the orientation of a vortex line. The plane repre- 
sents a section of the three-dimensional parameter space. The orientation is 
defined by the vorticity of the vortex which is the section of the vortex line. 
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Figure 3: An example for the neighborhood of a gap-closing point with 
det V > 0. The motion of the vortex line implies that the change of the 
Hall conductivity agrees with the change of the sign of the Dirac mass. 
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